Spin-wave velocities, density of magnetic excitations, and NMR relaxation in 

ferro-pnictides 
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We perform an analysis of the experimentally known temperature dependence of the staggered 
magnetization in the antiferromagnetic phase. This analysis allows us to put an upper limit on 
the unknown value of the spin wave velocity along the stripes of equal spin direction (spin stripes). 
The velocity is about ten times smaller than the velocity perpendicular to the spin stripes. The 
strongly anisotropic spin-wave dispersion implies a high density of low energy magnetic excitations. 
We demonstrate that this high density strongly enhances the ^^As NMR spin-lattice relaxation via 
the Raman scattering of magnons. We derive the polarization dependence of this relaxation channel 
and find very good agreement with experimental data. The high density of low energy magnetic 
excitations deduced from our phenomenological analysis supports the scenario that ferro-pnictides 
are close to a quantum phase transition. 
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I. INTRODUCTION 

One of the most important and widely discussed is- 
sues in the physics of iron pnictidc superconductors is 
whether these materials are strongly or weakly corre- 
lated. A closely related issue is the origin of magnetism in 
their parent compounds. In a slightly simplistic way one 
can formulate the problem in the following way. Does the 
magnetism arise from itinerant electrons or is it due to 
localized electrons? For a recent review, see for instance 
Ref. Q. The situation is different from cupratcs where 
parent compounds are clearly Mott insulators and hence 
there is no ambiguity about the origin of magnetism. 

In the present work, we do not address the issue of 
strong or weak correlations directly. In a phenomeno- 
logical way we analyse available experimental data on 
low temperature magnetic properties and determine the 
previously unknown spin-wave velocity along the spin 
stripes, by which we refer to chains of spins running in 
h direction in which the spins point all in the same di- 
rection (cf. Fig.[T]). This spin wave velocity turns out to 
be very small. It is by an order of magnitude smaller 
than the velocity perpendicular to the spin stripes. The 
knowledge of the velocity is very important itself because 
it predicts the outcome of future inelastic neutron scat- 
tering measurements. In addition, the knowledge sheds 
light on the issue of strong or weak correlations. The low 
velocity implies a high density of magnetic excitations 
and the high density strongly supports the strong corre- 
lation scenario based on the vicinity to a quantum criti- 
cal point^i^ii^. Due to the high spectral density magnons 
must contribute significantly to the NMR relaxation rate 
at a temperature above the spin-wave gap. We consider 



this mechanism for ^^As NMR spin-lattice relaxation and 
find a very good agreement between theoretical results 
and experimental data. 

The magnetic long range order is firmly established in 
the parent compounds LaFcAsO and Sr(Ba,Ca)Fe2As2 
by neutron scattering^i^i^, muon spin resonance, and 
Moessbauer spectroscopyiSiii. The neutron scattering 
reveals a columnar antiferromagnetic ordering with a 
staggered magnetic moment of (0.3 — 0.4)^b in LaFeAsO 
and (0.8 - 0.9)^b in Sr(Ba,Ca)Fe2As2. All the com- 
pounds are layered systems consisting of Fe-As planes. 
For simplicity, we consider only the tetragonal lattice 
which is formed by the Fe ions ignoring a small or- 
thorhombic and even monoclinic structural distortion. In 
Fig. [1] we show schematically the Fe-As plane and the 
spin ordering at the Fe sites. Along the a axis the spin 
directions alternate whereas they are the same along the 
b axis. Spins also alternate along the c axis which is or- 
thogonal to the plane. In our study, the lattice spacings, 
i.e., the distances between Fe ions, are ~ <?6 ~ 2.79A 
and (7c ~ 6.15A. We choose units such that all lattice 
spacings equal unity, ^ 1, 1, ^ 1- Note that 

the arsenic ions are shifted out of plane by 5c ~ ±1.35 A 
in a checkerboard pattern shown in Fig. [TJ 

The spin wave velocities along the a and the c axis 
as well as the spin wave gap at zero temperature have 
been measured by neutron scattering for SrFe2As2^ and 
for BaFe2As2^, 



A(T 



Va « 205meV , 
Vc ~ 45meV , 
: 0) w 6.5mcV . 



(1) 
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But the spin- wave velocity along the b axis, i.e., along 
the spin stripes, has not yet been measured to our knowl- 
edge. Note that the unit cell lattice spacings for the 122 
compounds SrFe2As2 and BaFe2As2 are twice larger than 
the corresponding values of g, a = 2ga, b = 2gh, c = 2gc- 
The standard crystallographic convention is to set a — > 1, 
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FIG. 1: The Fe-As plane. Fe ions are shown by open circles 
and As ions are shown by filled circles. The Fe spins are shown 
by arrows. The Fe ions lie exactly in the plane while As ions 
are out of plane by Sc « ±1.35 A in a checkerboard pattern. 
The pattern is shown by different fillings of the symbols for 
the As ions. 



— > 1, c ^ 1. Thus the values of the spin wave velocities 
in these standard units are twice larger than the values 
in our units. 

The Neel temperature for these compounds is Tjv = 
200 — 220K. The temperature dependence of the normal- 
ized intensity of elastic neutron scattering, /(T)//(0), 
and the temperature dependence of the normalized spin- 
wave gap, A(T)/A(0) have been measured in Ref. 0. 
These experimental results are shown in Fig. O 

We are not aware of a direct measurement of the spin 
magnetic susceptibility of SrFe2As2. However, the data 
for LaFeAsOi^iH and BaFe2As2ii show that the spin sus- 
ceptibility in the AF ordered phase averaged over direc- 
tions is about 



Xs 



1 X 10 *emu/(mol Fe) 



(2) 



We will use this value for SrFe2As2 having in mind that 
it might be by a factor « 1.5 — 2 off. 



II. EFFECTIVE ACTION FOR MAGNETIC 
EXCITATIONS 

To describe spin waves we use an effective model, the 
nonlinear cr-model with the following Lagrangian 

L = 

(3) 



where 



df_t = {dt, ivada, ivbdb, ivcdc) 



(4) 



and Aj is the bare spin-wave gap. The standard con- 
straint = 1 is imposed. Hereafter we set fc^ = h = 1 
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FIG. 2: (color online) Temperature dependence of the in- 
tensity of elastic neutron scattering at the antiferromagnetic 
superlattice reflection and of the spin- wave gap in SrFe2As2. 
Points with error bars show experimental data from Ref. 0. 
Blue circles show the normalized neutron scattering intensity, 
/(r)/J(0), and red diamonds show the normalized spin- wave 
gap, A(T)/A(0). The curves show theoretical results for the 
normalized neutron scattering intensity for various values of 
the spin-wave velocity Vb along the spin stripes. Solid black 
curves correspond to the first scenario for the spin-wave gap, 
Eq. (|25p . and dashed red curves correspond to the second 
scenario, Eq. (|26|) . The theory is justified only where the 
deviation of I(T)/I{Q) from unity is small. 



for simplicity. Note that generally the bare spin-wave 
gap Ah can depend on temperature if the temperature 
dependence arises from physics different from spin waves, 
say from phonons. We stress that this phenomenological 
description of low energy magnetic excitations is valid 
independently of the specific mechanism for magnetism. 
The description is equally valid for magnetism caused 
by itinerant electrons and for magnetism caused by lo- 
calized electrons. We will use ((3)) below the Ncel tem- 
perature T/v. The field theory ^ is the only possible 
effective theory that describes spin waves with disper- 



sion ujq = \/v\q^ + v'lq^ + t;2q2 ^ Therefore, the 



only important issue for the justification of ([3]) is that 
there are well-defined low-energy spin waves. This is di- 
rectly supported by experimenl^ii^. According to Ref. H 
the spin waves are well-defined up to w « 150meV. This 
is an important piece of information, but pragmatically, 
for purposes of the present work, we only need that spin 
waves exists with energies uj ^ Tn « 20meV. 

The spin-wave velocities Va and Vc are known from 
experiment, see Eq. ([T|). The susceptibility x± is related 
to the spin magnetic susceptibility ([2]) 



(5) 



where g is the gyromagnetic ratio, /^s is the Bohr mag- 
neton, Na is the Avogadro constant, and the factor 2/3 
comes from averaging over orientations. Lacking any 
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other information, we take the standard value of the gy- 
romagnetic ratio, g = 2. Eqs. ([5]) and ^ yield 



X± = 1.2 • lO^^l/meV. 



(6) 



Due to the uncertainty in the values for Xs and g one has 
to face an uncertainty in the value of x±- At worst, we 
estimate the uncertainty in the value of x± to be a factor 
1.5-2 relative to the value given in Eq. 

Assuming that the system is below the Neel tempera- 
ture, fi « (na,0,0), we represent the staggered magneti- 
zaton by 



n = {na,nb, n^) = (Vl - tt^, tt) 



(7) 



where the field tt is two-dimensional having only b- and 
c-components. The static component of the staggered 
magnetization reads in the first two leading orders 



(8) 



where (. . .) denotes the quantum Gibbs expectation 
value. 

Expanding Eq. ([3]) in powers of tt up to quartic terms 
(single loop corrections) we obtain the following La- 
grangian for the Tr-field 

L^^{{d,n)'+n\d,nr~Aln'} . (9) 

In order to derive the temperature dependent quadratic 
effective Lagrangian Lt one has to perform a decoupling 
in the quartic term in ^ 



(10) 



Here we have used integration by parts and the equation 
of motion in leading order 



(11) 



Hence the decoupling yields the following effective La- 
grangian 

= ^ {[1 + {T^'md.nf Al[l {7:')]n^} . (12) 



After rcscaling the field, ttji = -y/l + (tt^)??, wc find that 
the spin- wave gap is renormalized as 



A(T) = (l-(Tf2))A, 



(13) 



Following Fermi's Golden Rule, the scattering intensity 
I(T) is proportional to the square of the staggered mag- 
netization (ria)^ which is given by Eq. ([8]). Combining 
this fact with Eq. (fT3|) yields 



l-(7f^)+0((7r^)^) , (14) 



i{T) _ A(r) 



where is a bare scattering intensity without quantum 
fluctuations, for further discussion see below. In deriving 
this equation we assume formally that (7?^) ^ 1. Wc will 
discuss this point in more detail below. 

A standard calculation of the expectation value (tt^) 
leads to the following result 



1 



d^q 1 
X± J (27r)3 ujg 



(2n(^,) + 1) 



where 
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A2(r) 



is the spin-wave dispersion and 



n{ujq) 



1 



(15) 



(16) 



(17) 



is the Bose-Einstein distribution function. The unity in 
the factor (2nq + l) in (|15p is due to quantum fluctuations 
which lead to the quantum renormalization of the bare 
spin wave gap, Ab{T) Aqr{T), and to the concomitant 
renormalization of the bare scattering intensity, Ih 
Iqr (T) ■ The quantitative outcome of this renormalization 
depends on the high-energy cutoff 



h 



A 



qr 



d^q 1 
X± J (27r)3 ujg ' 



(18a) 



In our study, wc include the case that the quantum renor- 
malized quantities retain a temperature dependence from 
a temperature dependent Ab{T), whose dependence is 
induced from physical effects outside of the non-linear a- 
model, for instance from structural changes or phonons. 
The quantum renormalized quantities equal the physi- 
cal ones at zero temperature A(0) = Agr(O) and /(O) = 
lqr{0). Note that lb is temperature independent by defi- 
nition. 

On the present single loop level, the temperature ef- 
fects can be accounted for by 



I{T) A(T) 



Igr{T) Agr{T) 



1 - (t^'] 



the IT 



(19a) 



therm 



X± Jo 



p{uj,A(T))n{uj)duj (19b) 



instead of Eqs. p4)l and pS]) . Here we use the density 
of magnetic excitations p{uj, A) which reads in three di- 
mensions 

p{uj,A)^^^ wVt^2_ A2 e(cj-A) . (20) 

TT'^VaVbVc 

The above continuum expressions in three dimensions are 
only valid if Vb is not very small compared to temper- 
ature, VbTT ^ r, because for dominating temperature 
the boundaries of the Brillouin zone are felt which are 
not captured by the non-linear a model. In the opposite 



4 



limit. UbTT ^ T one should use the the two-dimensional 
density of 



1 



TTVaVc 



-UJ Q{uj - A) 



(21) 



We point out that for temperature independent bare 
gap Afe we have Iqr{T) = /(O) and A,^(T) = A(0). Then 
Eq. (|19ap already provides the result to be compared with 
experiment. Both normalized quantities, gap and inten- 
sity, should display the same temperature dependence in 
their deviation from unity. For this reason they are de- 
picted in the same plot in Fig. [2] We will discuss the 
very different behaviour of both experimental quantities 
below. 

Eq. (|19ap does not yet provide the ratio I{T)/I{Q) 
given by experiment if A^ is temperature dependent. To 
obtain full knowledge about /(T)//(0) we have to ac- 
count for the influence of the infrared cutoff, i.e., the 
gap, on the quantum renormalization. We find 



= 1 - 



(22a) 



V(0) 

1 f^^^pi.,A,AT))~pi.,A,M)^^,,^^ 



2x 



where we introduce a high-energy (UV) cutoff A to en- 
sure convergence. A realistic estimate is A = 200meV— 
Note that the difference occurring in (|22[) depends only 
weakly, i.e., logarithmically on the precise value of A. In 
two dimensions, there is even no dependence on the UV 
cutoff at all. 

The final result is obtained by combining (|19ap and 
(1221) in 



IjT) _ I{T) V(T) _ . 
7(0) Igr{T) V(0) L 



;^') therm] Rqr{T) (23) 



where /<jr(0) = /(O) entered. 

A remark on the validity of the single loop approxima- 
tion is in order. Obviously, the theoretical expressions are 
only valid if the thermal renormalization remains small, 
i.e., SI/I{0) < 1, (5/ 7(0) - I{T). Therefore, one can- 
not rely on (|23p in the vicinity of the critical point. But 
we can rely on at T < 200K where, according to 
experimental data, (57/7(0) < 0.3. 



III. VALUE OF THE GYROMAGNETIC RATIO 
AND POSSIBLE ORBITAL DYNAMICS 

In the present work we use the standard value g = 
2 which is the spin gyromagnetic ratio. This scenario 
assumes that the magnetism in the system is entirely 
due to spins. But the orbital physics of iron is certainly 
more complex and spin-orbit coupling plays an important 
rolei^. So values of the g-factor or g-tensor different from 
a scalar value of 2 are well possible. 



The phenomenological nonlinear cr-model description 
is valid independent of the origin of magnetism. There- 
fore our conclusions in the present paper do not depend 
on the extent that orbital and/or charge degrees of free- 
dom play a major role. Only the precise numerical es- 
timates depend on the numerical value of g. So for the 
present paper a better knowledge of g will affect only the 
numerical estimates, not the scenario. But for micro- 
scopic considerations, for instance the issue which value 
of spin is most appropriate, the local orbital physics is 
of fundamental importance and measurements of g can 
shed light on this issue. 

We are not aware of data on the value of g. Hence we 
would like to point out that the large value of the spin- 
wave gap, A(0) = 6.5mcV, gives a unique opportunity to 
measure g by inelastic neutron scattering. If a magnetic 
field B is applied that is directed along the a axis the 
spin wave excitations will be shifted in energy according 
to their projection of the angular momentum ±1 along 
a. Therefore, the magnetic field will split the spin wave 
gap in two 



A_ = 6.5meV — .g^sS 
A+ = 6.5uicV + gfigB 



(24a) 
(24b) 



For B = 15Tesla and for g = 2 the splitting takes the 
value 2gpBB w 3.5meV so that it should be easily ob- 
servable in neutron spectra. 



IV. ANALYSIS OF EXPERIMENTAL DATA 

For explicit calculations we need to specify the temper- 
ature dependence of the quantum renormalized gap Ag^- 
Since we do not have experimental knowledge about this 
quantity we study two scenarios which correspond to op- 
posite limits. It will turn out that our conclusions depend 
only weakly on which scenario is realized. 

The spin wave gap is caused by spin-orbit interaction 
in combination with the orthorhombic lattice deforma- 
tion. The deformation is practically temperature inde- 
pendent below Tn. Hence the scenario (i) assumes that 
the quantum renormalized spin wave gap A^^ is temper- 
ature independent. Then the observed spin gap acquires 
its temperature dependence A(T) solely from Eq. (|19ap . 
The zero temperature value is fixed to 



A(0) = 6.5meV 



(25) 



and the temperature dependence of the physical gap 
A(T) is determined by the self-consistent solution of Eqs. 
([J), dll), (fT9a| . (pb]) . and ^ or However, accord- 
ing to Eq. (|19ap this scenario implies the identical tem- 
perature dependence of the normalized neutron intensity 
and of the normalized spin- wave gap. This consequence 
is not supported by experiment. According to the data 
from Ref. shown in Fig. [5] the dependencies are signif- 
icantly different. But it cannot be excluded that experi- 
mental difficulties, for instance the influence of the charge 
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degrees of freedom, prevent the reliable measurement of 
the spin gap at finite temperature. So it is instructive to 
consider scenario (i) as one limiting case. 

In scenario (ii) we assume that the quantum renormal- 
ized gap A^^ is temperature dependent in precisely such 
a way that the experimentally observed A(T) shown in 
Fig. [2] is induced. The temperature dependence Abii{T) 
may result from the influence of low-energy phonons. In 
this case we fit the experimental datai for A(T) by the 
linear function 

A(T) = 6.5meV - 0.020TmeV/K . (26) 

Then we employ Eqs. H]), dH), (fTga|) . (pb)) . (^0]) or (^1]) to 
calculate the neutron scattering intensity resulting from 
the phenomenological spin- wave gap (j26p . 

The theoretical results obtained in the two scenarios 
for the values of the spin-wave velocity along the spin 
stripes Vb = lOOmeV, 50meV, 30meV, 20meV, lOmeV, 
and OmcV arc displayed in Fig.O The solid black curves 
correspond to scenario (i), i.e., Eq. The dashed red 

curves correspond to scenario (ii), i.e., Eq. p6|) . 

We emphasize that both scenarios for the temperature 
dependence of the spin- wave gap yield almost coinciding 
curves for the neutron scattering intensity. This can be 
attributed to the fact that by construction both scenarios 
are equal at T = so that the difference between them 
can only be discerned at sufficiently large temperature. 
But if T ~ ^{T) the precise value of the spin wave gap 
does not matter anymore. 

We already pointed out that the comparison between 
theory and experimental data makes sense only below 
180-200K where SI/I{Q) is sufficiently small and the the- 
ory is quantitatively reliable. The curves with Vb = 
lOOmeV and Vb = 50meV clearly disagree with experi- 
ment, while the curves with Vb = 30meV and Vb = lOmeV 
constitute upper and lower bounds to the experimental 
data. The curve for Vb ~ 20mcV is in good agreement 
with the data in the range of its validity {51 /I <C 1). 
The 3D formula (|20|) is at the verge of its validity for 
Vb = lOmeV since Vbir ~ T. Therefore, in Fig. [2] we also 
include curves for Vb = Q which are obtained using the 
2D formula (PT|) . We also studied the 3D 2D crossover 
empirically replacing v^ql by sin^(gf,/2) in Eq. ([16]). 
We do not show the corresponding curves because they 
completely confirm the results shown in Fig. [21 

Thus, our conclusion is that the value of the spin- wave 
velocity along spin stripes, Vb is in the range 

Vb~lQ- 30meV, (27) 

that means it is at least ten times smaller than Va ~ 
205meV. We stress that in essence this conclusion is 
based on the density of magnetic excitations (|20p . The 
number of excitations at low energies govern the ther- 
mally induced reduction of the staggered magnetic mo- 
ment. If Vb is large, the density is low implying a weak 
temperature dependence of the elastic neutron scatter- 
ing intensity. If Vb is small, the density is large implying 



a strong temperature dependence of the elastic neutron 
scattering intensity. The value of Vb must be sufficiently 
low to produce the experimentally found temperature de- 
pendence of the scattering intensity. 

The very low value of Vb that follows from the ex- 
perimental result for I{T) implies that magnetic fluc- 
tuations play a prominent role in the ferro-pnictides. 
This supports the previously proposed scenario that the 
ferro-pnictides are systems close to a quantum phase 
transition^. 



V. NMR SPIN-LATTICE RELAXATION FOR 
^^AS. 

As an additional testbed for our scenario of strong 
magnetic fluctuations at low energies we study the NMR 
relaxation rate T^^. This relaxation is due to inelas- 
tic Raman type scattering of thermally excited magnons 
from nuclear spin, see Fig. [3| Clearly, this mechanism 




FIG. 3; Magnon Raman scattering on nuclear spin due to 
hyperfine interaction. Solid lines denote magnons, the dou- 
ble lines denote a nuclear spin. The dashed line denotes the 
hyperfine interaction. 

is most important at temperatures above the spin-wave 
gap. For concreteness, we consider the relaxation of the 
nuclear spin of ^^As. The As ion is positioned at the top 
of a pyramid with four Fc ions at its base, see Fig.[T| The 
hyperfine interaction of the ^^As nuclear spin / with the 
electronic spin S on the adjacent Fc ion is of the following 
formic 

H = B {I- S) + C {!■ N){S ■ N) , (28) 

where is a unit vector directed from As to Fe, see 
Fig. [Tj The B-term is due to s-wave transferred hyper- 
fine interaction and the C-tcrm is due to p-wave trans- 
ferred hyperfine interaction. It is knowni^ that the ef- 
fective static hyperfine magnetic field at As takes the 
value i?cff ~ 1.5Tesla and is directed along the c axis. 
The average electronic magnetic moment at the Fe site is 
« Q.S^B-^i^. Assuming that g ~ 2 this implies an average 
static spin component of (S) ~ 0.4. Since the electronic 
spins are arranged in a staggered pattern, see Fig. [TJ only 
the C-term in ([^S]) contributes to the effective static field. 
Hence, 

;i„Boff = 4 X 0.4 X 0.33 X C X / , (29) 
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where / = 3/2 is the nuclear spin and ^„ = 1.86^Ar is 
the magnetic moment of the ^^As nucleus. The factor 
4 in (|29p stems from the four neighboring Fc ions and 
0.33 = NaNc is a product of components of the unit 
vector N in the geometry of the As-Fe^ pyramid. From 
(I29ll we find 



C w 1.1 



10^'^mcV 



27MHz 



(30) 



Alternatively, the rescaling from the known values of 
transferred hyperfine constants in cuprates yields the fol- 
lowing estimates 



B 
C 



10 
10 



"^mcV 
""'mcV 



250MHz 
25MHz . 



(31a) 
(31b) 



While the estimate for C agrees very well with (|30p de- 
duced from the experimental dat a^^'^^i^^ , the estimate 
(|3T|) for B is about five times larger than the one mea- 
sured in Ref. It is worthwhile noting that there is a 
comment in Ref. [l^ that they might underestimate the 
value of B. In the present work, we will rely on the value 
of C given in (|30|) and on the estimate for B given in ([3T|) . 
It is very natural that B ^ C because B is due to the 
s-wave and C is due to the p-wave hyperfine interaction. 

The static components of the electron spins are polar- 
ized along the a-axis, see Fig. [TJ Hence the spin wave 
excitations are polarized along the b- and the c-axis. To 
describe the magnon Raman process shown in Fig. [3] we 
need only the part of (|28p that is bilinear in spin wave cre- 
ation and annihilation operators, i.e., bilinear in tt in the 
language of the non-linear a model. This implies that we 
only need to keep the terms in (pS)) that are proportional 
to the a-component of the electron spin 5* 



(32) 



where the summation goes over four nearest Fe sites. 
Finally, in the notation of the cr-model, S / Scs ft = 
{Ua, rib, Tic) = (Vl ^ TT^, tt), this leads to 

H ^ 5(r)5'eff(2 B la da + 4 0.33C /,)n, 

^ 6iv)SMB la da + 0.66C /Jtt^ , (33) 

where constant terms are omitted in passing to the last 
line. The gradient da along the a-axis in the B-term 
appears because the magnetization in this direction is 
staggered. We remind the reader that in our notations 
both S{r) and d are dimensionlcss. In Eq. ((33l) we have 
introduced the effective spin ^cff. The first naive impres- 
sion is that Scs = (S) « 0.4. This would imply that 
the magnon Raman operator is renormalized by quan- 
tum fluctuations exactly like the staggered magnetiza- 
tion. However, we have checked by an explicit single loop 
calculation that the Raman operator is not renormalized 
while the staggered magnetization is certainly reduced in 
the single loop approximation. So the naive expectation 
is wrong. For numerical estimates we will use 



5eff ■ 



1 



(34) 



It is clear from the kinematic structure of Eq. pS)) that 
the _B-term contributes to the spin-lattice relaxation only 
if the initial nuclear spin is directed perpendicular to the 
a-axis. The C-term contributes to the relaxation only if 
the initial nuclear spin is directed perpendicular to the 
c-axis. The Raman relaxation rate due to the _B-tcrm is 
given by Fermi's Golden Rule 



Wb 



2n 



y2 



(35) 



+ nq)5{ujk -tOq- Wnmr) ■ 



The factor 27r = 2 x ^ x 27r is the factor 27r from Fermi's 
Golden Rule multiplied by the number of magnon po- 
larizations, 2, and multiplied by i resulting the from 
Clebsch-Gordon coefficients related to the nuclear spin 
/ = 3/2. The factor (xlt^fet^,)^^ is due to the normal- 
ization of the 7? field, the factor {ka — qaY due to the 
gradient 9a, and Uq is given by p7|) . Since the NMR 
frequency is very small, wnmr <C ujq , the expression ([35]) 
can be transformed to 



Wb - 



[S,i,Bf 



12TT^{xi_VaVbVcf vl smh^{uj/2T) 



(^2 _ A2)2 



■du , (36) 



where A is the spin wave gap. The two factors (oj^ — A^) 
in the numerator of the integrand stem from the density- 
of-states in three dimensions and from the matrix element 
(fca - qaf- 

Similarly, the Raman relaxation rate due to the C-term 
reads 



Wc 



[5cff0.66C] 



(u;2 - A2) 



^-^^{xi-VaVbVcY J A smii{uj/2T) 



dcu , (37) 



where there is one factor (w^ — A^) less in the numera- 
tor because there is no particular momentum dependent 
matrix element. 

It is clear that Eqs. ([55]) and ([57]) are not justified in 
the vicinity of the Neel temperature. Obviously, they 
arc not valid at T > T/v either. So we use them only 
below Tjy. Both Wb and Wc are very steep functions 
of temperature. Plots of Wb and Wc calculated with 
parameters given by ([T|) , © , (HJ) , and Vb = 20meV 
are presented in Fig. U) The decay rate Wc is much 
smaller than Wb, Wb » Wc, and hence 1/Ti w Wb so 
that the magnon Raman relaxation is mainly due to the 
s-wave transferred hyperfine interaction. The estimate 
for the value of relaxation rate presented in Fig. 2] agrees 
with the available datai^iiiii^. 

One can fit the data by fine tuning of B and/or Vh 
and/or x±- But this is not our aim here. Note, how- 
ever, that the NMR relaxation rate (|36p is proportional 
to the second power of the density of magnetic excita- 
tions . So it is very sensitive to the value of p in (^0]) 
and hence to the value of Vb- Hence the fact that our 
previous set of parameters yields the correct magnitude 
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FIG. 4; NMR relaxation rates due to magnon scattering. 
The solid curve gives contribution of the s-wave transferred 
hyperfine coupling, see Eq. p6)) . The dashed curve gives con- 
tribution of the p-wave transferred hyperfine coupling, see 
Eq. (|37l). The parameters are given by ([TJ, ((6}, (|3l]|, ((SiJ and 
U6 = 20meV. 

of the NMR relaxation rate strongly supports our claim 
that Vh is small. 

Since the relaxation is dominated by Wb wc pre- 
dict a significant polarization dependence of the relax- 
ation. The corresponding Hamiltonian contains only the 
a-component of the nuclear spin, see Eqs. ([5^ and ([55)1 . 
Therefore, this mechanism does not contribute to relax- 
ation if the ^^As nuclear spin is polarized along the a- 
axis. The mechanism contributes equally for polariza- 
tions along the h- and the c-axis. For a twinned sample, 
where domains with swapped a and h axes are of equal 
weight, this argument implies that the relaxation for the 
c-polarization of nuclear spin is twice faster than the re- 
laxation for an in-plane polarization of the nuclear spin. 

For temperatures below the spin-wave gap the magnon 
Raman relaxation is essentially switched off. In other 
words, the mechanism related to collective magnetic 
modes is not active. But there is also a diffuse magnetic 
relaxation stemming from the charge degrees of freedom 
because the system is not an insulator. It is natural to 
assume that this diffuse relaxation scales linearly with 
temperature as it does in normal Fermi liquids. This 
low-temperature behavior of the relaxation was observed 
in Refs. [l^ and [H. The charge driven relaxation is 
certainly also active for T > A. But in this regime its 
contribution to the NMR relaxation rate is relatively 
small with the contribution from collective magnetic 
modes prevailing. 



VI. CONCLUSIONS 

We studied the parent ferro-pnictidcs below their Neel 
temperature. Based on the nonlinear cr-model, we consid- 



ered the available experimental data on the temperature 
dependence of the staggered magnetization phenomeno- 
logically. We found that one needs a high density of 
magnetic excitations to explain the relatively strong tem- 
perature dependence of the magnetization. This implies 
that the spin wave velocity along the spin stripes is very 
small. The values for this velocity estimated from the 
analysis are « 10 — 30meV. For comparison, the in- 
plane velocity perpendicular to stripes takes the value 
Va « 205meV. 

We also analyzed the NMR spin-lattice relaxation rate 
for ""^As. Due to their high spectral density the magnons 
dominate the relaxation rate at temperatures above the 
spin-wave gap. Our estimates for the relaxation rate 
based on the density found from the neutron scatter- 
ing data agree very well with direct NMR measurements. 
This is an independent confirmation of the high spectral 
density of magnetic excitations. 

So both the temperature dependent magnetization as 
well as the NMR relaxation rate confirm strong mag- 
netic fluctuations at low energies. Thus the present 
phenomenological analysis corroborates the scenario that 
the ferro-pnictides constitute systems close to a quan- 
tum phase transition triggered by frustrated magnetic 
couplings^. Hence a strongly correlated picture of the 
ferro-pnictides is favored. 

Very recent inelastic neutron scattering dat a^^i^" indi- 
cate the ratio of spin wave velocities Wa/fb ~ 2. This 
is not consistent with our conclusion Va/vb ~ 10. Our 
analysis of the temperature dependence of the staggered 
magnetization and especially of the NMR relaxation rate 
is in essence based only on the spin- wave dispersion (jl6p . 
Only the dispersion determines the density of excited 
magnons at a given temperature, and only the density 
determines the NMR relaxation rate. With the spin- wave 
velocity Vb taken fro m^^i^° one obtains the relaxation rate 
about 20-30 smaller than the experimental one. 

How can the above discrepancy be explained? The 
data— is taken on twinned samples because only be- 
low the structural transition temperature the orthorhom- 
bicity occurs. Superposing dispersions with prominent 
ridges such as the ones in Fig. 3a in Ref. can lead to 
responses similar to the ones in Ref— for moderate 
and high energies. The time-of-flight technique used in 
both experimental probe a^^'^" is certainly best suited for 
investigating the moderate and higher energies. 

If the careful study of the influence of twinning does 
not solve the discrepancy our analysis indicates the ex- 
istence of some low-energy (k, 10 — 20meV) magnetic 
degrees of freedom which have so far not been taken into 
account. These degrees of freedom must contribute to 
the NMR relaxation and they must be difficult to detect 
by neutron scattering. 

From our results and the above discussion we conclude 
that further experiments focusing on low lying magnetic 
modes are called for to resolve this crucial issue. It would 
be highly desirable if low-temperature detwinned samples 
could be generated. 
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